Density profiles of dark matter halos with 
anisotropic velocity tensors 



N. Hiotelis ^ 

1st Experimental Lyceum of Athens, Ipitou 15, Plaka, 10557, Athens, Greece, 
E-mail: hiotelis@avra.ipta.demokritos.gr 



Abstract 

We present density profiles, that are solutions of the spherical Jeans equation, de- 
rived under the following two assumptions: (i) the coarse grained phase-density 
follows a power-law of radius, p/cr"^ oc r~", and (ii) the velocity anisotropy param- 
eter is given by the relation Paif) = Pi + 2^2 i^^r/r )^ where P2 are parameters 
and equals twice the virial radius, r^j^, of the system. These assumptions are well 
motivated by the results of N-body simulations. Density profiles have increasing 
logarithmic slopes 7, defined by 7 = — The values of 7 at r = \^~'^'^rmr, a 
distance where the systems could be resolved by large N-body simulations, lie in the 
range 1. — 1.6. These inner values of 7 increase for increasing (3i and for increasing 
concentration of the system. On the other hand, slopes at r = r^r he in the range 
2.42 — 3.82. A model density profile that fits well the results at radial distances 
between 10~^r„jr and r^r and connects kinematic and structural characteristics of 
spherical systems is described. 
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1 Introduction 



Density profiles of galactic halos are well fitted by two-power law models. Nu- 
merical studies (Quinn et al. [1986]; Frenk et al. [1988]; Dubinski & Garlberg 
[1991]; Crone et al. [1994]; Navarro et al. [1997], NFW; Cole & Lacey [1996]; 
Huss et al. [1999]; Fukushige & Makino [1997]; Moore et al. [1998], MGQSL; 
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Jing & Suto [2000], JS) showed that the profile of relaxed halos steepens mono- 
tonically with radius. The logarithmic slope 7 = — is less than 2 near the 
center and larger than 2 near the virial radius of the system. The value of 
7 near the center of the halo is not yet known. Navarro et al. [1997] claimed 
7 = 1 while Kravtsov et al. [1998] initially claimed 7 ~ 0.7 but in their revised 
conchisions (Klypin et al. [2001]) they argue that the inner slope varies from 
1 to 1.5. Moore et al. [1998] found a slope 7 = 1.5 at the inner regions of their 
N-body systems. 

N-body simulations give also valuable information about the kinematic state 
of the relaxed structures. Taylor & Navarro [2001], (TN) showed that in a 
variety of cosmologies the velocity dispersion and the density p at distance 
r are connected by a scale free relation p/a^ oc r~" with a = 1.875. They also 
examined density profiles that resulted from solving the spherical Jeans equa- 
tion under the above condition assuming an isotropic velocity distribution. 
Final structures of N-body simulations are radially anisotropic. The radial 
anisotropy increases outwards taking its maximum value (that is in the range 
0.3 — 0.8) at about twice the virial radius of the the system. This behavior 
seems to be universal (Cole & Lacey [1996], Cohn et al. [2000], Carlberg et al. 
[1997]). 

In this paper wc present solutions of spherical Jeans equation using the above 
scale free relation given in TN and a model for the velocity anisotropy pa- 
rameter, in order to construct density profiles for radially anisotropic models. 
The method for the solution of Jeans equation in described in Sect. 2 and the 
results are presented in Sect. 3. The conclusions are discussed in Sect. 4. 



2 The solution of Jeans equation 



Under the assumption of spherical symmetry and no rotation the Jeans equa- 
tion is: 
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where is the velocity anisotropy parameter given by the relation 
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2^2 (^r) ('^h^re (7j and are the tangential and the radial dispersion of 
velocities respectively) and $ is the potential. Multiplying both sides of (1) 
with and differentiating with respect to r we have 
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Under the assumption 

p{T)/a\r) = (p(ro)/a3(ro))(r/ro)-^ 



(3) 



and after setting Pq = p(ro),o"o = <^{f^o),x = ln(r/ro) and y = ln(p/po) Jeans 
equation is written in the form of the foUowing set of differential equations: 



y2 
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with yi{x) = y{x) and 
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where k = ^ Note that cr in (3) is the total dispersion of velocities. 

TN showed that the isotropic case = 0) admits a solution yi{x) = —(3x, 
(power law in r — p space), with f3 = 6 — 2a and k = 2{j3 — 1)(3 — /9)/3. Ad- 
ditionally they used a = 1.875, a value resulting from the final structures of 
N-body simulations, to show that for values of k larger than the one that cor- 
responds to the above power-law solution, the density profiles become complex 
and that for values of k larger than a critical value kcrit, the density profiles 
become unrealistic since densities are vanished at some finite radius near the 
center of the system. The density profile that corresponds to kcru scales as 
j,-o.75 jj-^ ^Yie inner region of the system, close to the NFW profile that scales 
as r~^. 

Since the anisotropy of velocities is always present in the results of N-body 
simulations, it should be interesting to study its effects on the solutions of 
Jeans equation. N-body simulations for a variety of cosmologies show an al- 
most universal variation of /3a. Starting from the center of the system, 
increases outwards taking its maximum radius at about twice the virial radius 
of the system. The maximum value is in the range 0.3-0.8 (Cole & Lacey 
[1996], Cohn et al. [2000], Carlberg et al. [1997]). Thus we model /3„ as 



where r is the radial distance. The maximum value of Pa is Pmax — Pi + P2 
and it is achieved at r = r^, = 2r^ir. The initial conditions required for the 
solution of the above set of differential equations are set in a similar way as 
in TN: At a radius Tq = 1 we assume a density po = 1 and 7(ro = 1) = 2.25. 
Thus we take 2/1 (0) ~ and y'i{0) — —2.25. Furthermore, the value of 1/2(0) is 



3 



calculated by solving equation (5) for y2 leading to the relation: 
1 



Z/2(0) = 



3 - /5a(0) 



Y - Y + mm - Warn-' + 2/3a(0) 



(7) 



Obviously ^2(0) ~ Then we integrate inwards and outwards to calculate 
density profiles. 



3 Density profiles 



The calculation of density profiles is organized as follows: For given values 
of ryir,(3i,P2 we find the solution of Jeans equation for the maximum value 
of k, named kcrit, for which the density is a decreasing function of the radial 
distance. An example of the behavior of the solutions in the inner region of 
the system is presented in Fig. 1. The curves shown in the left panel are 
solutions for j3i = 0,^2 = 0.8 and rmr/ro = 10. The upper line corresponds to 
k — 0.84, the middle one is the critical solution corresponding to A; = kcrit — 
0.85114 and the lower one corresponds to k — 0.86. The middle panel shows 
solutions for r^ir/ro — 5, (3i — and /32 = 0.3. The upper line corresponds to 
k = 0.85, the middle one is the critical solution corresponding to k = 0.86146 
and the lower one to k = 0.87. The right panel contains the solutions for 
Pi = 0.15, = 0.35 and ryir/rg = 10. The critical solution (middle curve) is 
obtained for k — kcrit = 0.80031 while the upper and lower curves are solutions 
for k — 0.79 and k — 0.81 respectively. It is shown in this Fig. that for values 
of k smaller than kcrit the slope of the density profile changes significantly at 
the inner region of the system where it becomes very steep, while for k > kcrit 
density profile becomes unrealistic since it vanishes at a finite distance from 
the center. Although there is no obvious physical reason to exclude solutions 
with k < kcrit we note that such a behavior of the density profile has not 
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Fig. 1. Density profiles derived from the solution of Jeans equation for fixed values 
of Pi,P2 and c and varying k. For values of k larger than a kcrit densities vanish at 
a finite radius near the center while for k < kcrit the slopes increase inwards. 
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Table 1 

The parameters and the resulting values of kcrit for Cases 1-21 



Case No 


c 


Pi 


P2 


kcrit 


Case 


No 


c 


/3i 


P2 


kcrit 


Case 1 


10. 


0.00 


0.80 


0.85114 


Case 


10 


10. 


0.15 


0.65 


0.78410 


Case 2 


5. 


0.00 


0.80 


0.80947 


Case 


11 


5. 


0.15 


0.65 


0.74800 


Case 3 


2.5 


0.00 


0.80 


0.72034 


Case 


12 


2.5 


0.15 


0.65 


0.67110 


Case 4 


10. 


0.00 


0.50 


0.86638 


Case 


13 


10. 


0.15 


0.35 


0.80031 


Case 5 


5. 


0.00 


0.50 


0.84107 


Case 


14 


5. 


0.15 


0.35 


0.78150 


Case 6 


2.5 


0.00 


0.50 


0.78992 


Case 


15 


2.5 


0.15 


0.35 


0.74380 


Case 7 


10. 


0.00 


0.30 


0.87642 


Case 


16 


10. 


0.15 


0.15 


0.81100 


Case 8 


5. 


0.00 


0.30 


0.86146 


Case 


17 


5. 


0.15 


0.15 


0.80310 


Case 9 


2.5 


0.00 


0.30 


0.83184 


Case 


18 


2.5 


0.15 


0.15 


0.78750 












Case 


19 


10. 


0.30 


0.50 


0.71704 












Case 


20 


5. 


0.30 


0.50 


0.68700 












Case 


21 


2.5 


0.30 


0.50 


0.62330 



been reported elsewhere and so we restrict ourselves to the study of critical 
solutions alone. 

We studied 21 different cases. The values of the parameters for every case 
as well as the resulting values of kcrit are given in the Table 1. A convenient 
concentration parameter c is defined by the relation c = Vyir/ro, that is the 
ratio of the virial radius to the radius where the logarithmic slope of the 
density profile equals to 2.25. 

TN illustrated that the critical solution (k — kcrit) for the isotropic velocity 
case yielded a phase-space distribution function {dM/dlog{(T^/p)) that was the 
most peaked, or 'maximally mixed' (excluding cases that produced hollow-core 
density profiles). We show that the same holds in anisotropic cases. The left 
panel of Fig. 2 presents the phase space distribution function for the case 
with c = 10, 6i = and 62 = 0.8. The most peaked curve corresponds to 
k = kcrit = 0.85114 the middle one to A; = 0.84 and the lower one to A; = 0.82. 
The right panel presents the case with c = 10, &i = 0.15 and 62 = 0.35. The 
most peaked curve of the this panel corresponds to k = kcrit = 0.80031, the 
middle one to /c = 0.78, while the lower one to A; = 0.76. 

The resulting density profiles are fitted by curves of the form 
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Models with density cusps that have been proposed in the hterature belong 
to the above class of density profiles. For A = = 1 and u = 2 Eq.(8) 
corresponds to the NFW model. The model proposed by Hernquist ([1990]) 
has X — 1, pL — 1 and v — 3. MGQSL proposed the model with A = 1.5, pL — 1.5 
and 1/ — 1 while JS estimated X — 1.5, /i — 1 and v — 1.5 for galaxy-size halos. 
The fitting parameters pc, rg, X, /i and u of Eq.(8) are calculated finding the 
minimum of the sum 

NP 

S = Y.[logpj{r,)-logpf,t{n)f (9) 
1=1 



where pj is the density profile predicted by the above described solution of 
Jeans equation. The minimum of S is found using the unconstrained minimiz- 
ing subroutine ZXMWD of IMSL mathematical library. 

We used values of the density at radial distances in the range to 
l.Sr^jr- The logarithmic slope of the density profile at r — > is given by the 
value of A. For the first nine cases, that are all isotropic at the center = 0), 
values of A that result from the above fitting are in the range 0.71 — 0.98. For 
the next nine cases (Cases 10 to 18) that have j3i = 0.15 the values of A are 
in the range 0.82 — 1.02 while for the last three cases (Cases 19, 20 and 21), 
that are highly anisotropic at the center {Pi = 0.3) the values of A are in the 
range 1.06 — 1.22. We also note that in all cases the values of p are in the 
range 0.32 — 0.61 in contrast to all the proposed models that give p>l. 



The quality of fit is excellent and it allows a very good estimation of the log- 
arithmic slope of the density profile all the way from very small to vary large 
radii through the relation: 




Fig. 2. The phase-space distribution function for anisotropic cases. The most peaked 

curves correspond to the critical values of k for every case. The left panel presents 
systems with c = 10 6i = and 62 = 0.8 while the right one shows systems with 
c = 10 61 = 0.15 and 62 = 0.35. The systems of every panel have the same mass 
and energy. 
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7 = 



dln(p(r)) 



— \-\- jiv 




(10) 



dlnr 



1 + {r/r-sY 



The results for all cases are presented in Fig. 3 and Fig. 4. Figure 3 consists 
of six panels that show the density profiles as well as the logarithmic slopes 
of these profiles for the first nine cases. The two figures of the first column 
correspond to the Cases 1, 2 and 3 (sohd, dotted and dashed lines respectively). 
The second column shows the results for the Cases 4, 5 and 6 respectively while 
the third one contains the results for Cases 7, 8 and 9. Densities are normalized 
to the virial density of every case, that is the density of the system at its virial 
radius, while distances are normalized to the virial radius. In a similar way 
the rest 12 cases are presented in Fig. 4. The first column shows the Cases 10, 
11 and 12 (solid, dotted and dashed line respectively). Cases 13, 14 and 15 
are shown in the two figures of the second column. The third column shows 
the Cases 16, 17 and 18 and the last one corresponds to the three last Cases 
19, 20 and 21. In this way every column shows three cases that have the same 
values of /3i and (32 but different values of c. Additionally, irrespectively of the 
column, solid lines correspond to systems with c = 10, dotted lines to systems 
with c — h and dashed lines to systems with c—2.h 

The results lead to the following conclusions: 



Fig. 3. Density profiles and their logaritmic slopes for Cases 1 to 9. First column: 
Cases 1, 2 and 3 (solid, dotted and dashed line respectively). Second column: as in 
the first column but for Cases 4, 5, 6. Third column: as in the first column but for 
Cases 7, 8, 9. 




log( r / r„ J 



7 




-3 -2 -1 -3 -2 -1 

log(r/rJ 

Fig. 4. As in Fig.3 but for Cases 10-21. 



(1) The slope of the density profile at the central region is closely connected 
to the amount of the anisotropy of velocities there. It is larger for larger 
anisotropy. For example at r = 10~^-^r„jr the first nine Cases 1 to 9 that 
have Pi = 0.0 show central slopes near to unity [1 — 1.18]. Cases 10 to 18 
that correspond to f3i = 0.15 have slopes in the range [1.25 — 1.39] and 
the last three cases, 19, 20 and 21 that have highly anisotropic central 
regions, = 0.3), have slopes in the range [1.49 — 1.61]. 

(2) Structures resulting for the same values of Pi and P2 show inner slopes 
that depend on their concentrations. This dependence is not very clear 
at very small radii, (for example at r = lO^^'^r^jj,.) but it is obvious at 
larger radii as those that can be resolved by N-body simulations, (r ~ 
10~^'^r^ir)- Thus, if the values of Pi and P2 are universal, systems with 
larger c have larger slopes than those with smaller c at distances that 
are the same fraction of their virial radii. Note that the results of N-body 
simulations suggest that c is a decreasing function of the virial mass of the 
system (NFW, Bullock et al. [2001]). Thus, our results indicate that inner 
slopes are smaller for systems with larger masses. This dependence is in 
accordance with the results of N-body simulations. For example JS found 
that, at r = O.Olr^j^, the slope is ~ 1.5 for galaxy-size halos and ~ 1.1 
for galaxy cluster-size halos respectively. We note at this point that our 
concentration parameter c is related to the concentration parameter of 
NFW, cnfw = Tmr/Ts,NFWi by the relation cnfw = (3/5)c. The results 
of N-body simulations indicate that for galaxy-size halos c^fw ~ 17 (that 
corresponds to c = 10), while for galaxy cluster-size halos cj^fw ~ 7 that 
corresponds to c ~ 4.2. 

(3) Slopes at r = r^r show values in the range 2.42 — 3.82. For systems with 
the same values of Pi and P2 outer slope is an increasing function of c. On 
the other hand, a comparison between systems with the same values of 
c and Pi shows that the ones with smaller values of P2 have larger outer 
slopes. 

We summarize the above results by proposing a fitting formula (FF), that 
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holds for the range 10 "^Tmr to and connects the values of the exponents 
A, // and u of Eq. (8) with the values of /?i and /?2- That is: 

A = 1 + 2/3i, II = 0.75 - 0.3/3i, ^ ^ ^ - 2/3i - ^^^^ 

A* 



The two remaining fitting parameters pc and are found by the minimization 

procedure described above. The results are shown in Fig. 5. Solid lines are the 
solutions of Jeans equation while squares are the respective fits given by FF. 
The fit is more than satisfactory. 

Additionally, the resulting concentrations are compared to the input concen- 
trations by the following way: Consider r„ the radius where the slope equals 
to n. This r„ is related to the fitting parameters by the relation: 

^^={t^^]'^- (12) 
\X-\- pLU — nj 



The concentration of the system, resulting from the fitting procedure, is c/ = 
Tvir/f 2.25- The values of c/ are close enough to the input values of the con- 
centration c. For systems with c = 10 the values of c/ arc in the range of 
7.8 — 12.58, for systems with c = 5 the values of C/ lie in the range 4.28 — 5.33, 
while the density profiles of systems with c = 2.5 are fitted by density profiles 
with Cf in the range 2.02 — 2.8. 

In Fig. 6 a detailed comparison between various models is shown. For Cases 
1 and 5 the solutions of Jeans equation are shown as solid lines, the results 
of FF as dotted hne and the NFW profile as dashed lines. For cases 14 and 
18 sohds and dotted hues show the same profiles as in Cases 1 and 5 while 
dashed lines correspond to the MGQSL profile. It is clear that the NFW profile 
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Fig. 5. A comparison of all resulting density profiles with our model. Solid lines are 
the solutions of Jeans equation while the predictions of our fitting models are shown 
by squares. 
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is an excellent fit for the Case 5 and a very good fit for the Case 1. On the 
other hand MGQSL is an excellent fit to Case 14 and a good enough fit to 
Case 18. The FF gives very good fits for all four cases. It is characteristic 
that analytical models with different values of the parameters A, /x and v can 
fit equally well the same density profile. This indicates the large degree of 
degeneracy in these parameters that is already noticed by Klypin et al. [2001]. 
However a question natmally arises at this point. Can both popular models of 
NFW and MGQSL fit equally well the same density profile? The answer is no 
as it can be obtained by looking at Fig. 7. Sohd hues are the solutions of Jeans 
equation, dashed lines are the best NFW fit while dotted hues correspond to 
MGQSL fit. The results of FF are shown as small circles. For the centrally 
isotropic Case 8 the NFW profile is an excellent fit. For this Case the MGQSL 
profile does not approximate well the results. On the other hand the results of 
Case 10, that corresponds to a centrally anisotropic system, are fitted well by 
the MGQSL profile while the NFW is less steep near the central region and 
also differs at intermediate distances. The FF is a very satisfactory fit to both 
cases. 




Fig. 6. Comparison of resulting density profiles with other models. Solid lines: solu- 
tions of Jeans equation, dashed lines: NFW fit (for cases 1 and 5) and MGQSL fit 
(for cases 14 and 18), circles: FF fit. 
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Fig. 7. Comparison of resulting density profiles with other models. Solid lines: so- 
lutions of Jeans equation, dashed lines: NFW fit, dotted lines: MGQSL fit, circles: 
FF fit. 

4 Discussion 



This paper presents density profiles of spherical systems derived by the solu- 
tion of spherical Jeans equation. These profiles, derived under two assumptions 
described in previous sections, show a considerable agreement with models pro- 
posed in the literature. The results show a number of trends relating kinematic 
and structural characteristics of the systems. A formula that fits very well the 
results for r > lO^^Vyir and connects the values of the anisotropy velocity 
parameter to the form of the density profile is proposed. If all cases studied 
above do represent final states of real systems then the idea of a density pro- 
file with fixed values of \, fi and u able to fit all halos, by just adjusting the 
concentration, is not favored. However it should be interesting, large N-body 
simulations, that are the most powerful methods to deal with the formation of 
structures, to answer some of the questions that naturally arise by our results. 
Such questions relate to the validity of our two main assumptions, the univer- 
sality of the values of f3i and as well as the relation between the values of 
these two parameters and the mass of the system. Answering these questions 
should reduce the parameter space and help to construct better analytical 
models so as to improve our understanding about the formation of structures 
in the universe. 
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